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The shear (77) and bulk (f) viscous coefficients have been evaluated for a gluonic fluid. The elastic, 
99 ~* 99 an d the inelastic, number non-conserving, gg —¥ ggg processes have been considered as 
the dominant perturbative processes in evaluating the viscous co-efficients to entropy density (s) 
ratios. Recently the processes: gg — > ggg has been revisited and a correction to the widely used 
Gunion-Bertsch (GB) formula has been obtained. The 77 and ( have been evaluated for gluonic fluid 
with the formula recently derived. At large a s the value of n/s approaches its lower bound, ~ l/4n. 
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The nuclear collisions at Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider 
(LHC) energies are aimed at creating a phase where 
the properties of the matter is governed by the 
quarks and gluons [l|, such a phase, which is mainly 
composed of light quarks and gluons - is called quark 
gluon plasma (QGP). The weakly interacting picture 
of the QGP stems from the perception of the QCD 
asymptotic freedom at high temperatures and den- 
sities. However, the experimental data from RHIC, 
especially the measured elliptic flow [2[ indicate that 
the matter produced at Au+Au collisions exhibit 
properties which are more like strongly interacting 
liquid than a weakly interacting gas. The shear vis- 
cosity or the internal friction of the fluid symbolizes 
the ability to transfer momentum over a distance of 
~mean free path. Therefore, in a system where the 
constituents interact strongly the transfer of momen- 
tum is performed easily - resulting in lower values of 
i]. Consequently such a system may be character- 
ized by a small value of rj/s. The importance of 
viscosity also lies in the fact that it damps out the 
variation in the velocity and make the fluid flow lam- 
inar. A very small viscosity (large Reynold number) 
may make the flow turbulent. 

On the other hand the bulk viscosity exhibit the 
exchange of energy between the translational and in- 
ternal degrees of freedom. Although much emphasis 
has been given to the evaluation of the shear viscos- 
ity for a partonic system recently, the bulk viscosity 
is comparatively less discussed. Probably, because 
the bulk viscosity for a structureless point parti- 
cles vanishes both for relativistic and non-relativistic 
limits However, there are several reasons for 
which the bulk viscosity of a system formed in nu- 
clear collisions at ultra-relativistic energies may be 
non-zero 0]. The trace anomaly in QCD will give 
rise to non-zero £, which will indicate the deviation 
of the system from the conformal invariance, because 
the C is defined as the correlation of the trace of the 
energy momentum tensor through Kubo's formula. 
The £ for SU(3) gauge theory has been evaluated 
in lattice QCD and its value is found to be quite 



large around the temperature domain for partons to 
hadrons transition (T c )[5[. The divergence of £ may 
be treated as a signal of critical point as it diverges 
near this point [6|. However, this point has been 
challenged in Q- This indicate that both 77 and ( 
can be used effectively to characterize QGP. There- 
fore, in the present work we would like to estimate 
both the ratios, rj/s and £/s for a gluonic fluid by 
taking into account the perturbative QCD (pQCD) 
elastic process, gg — > gg and the inelastic, number 
non-conserving process, gg — > ggg |8j. While evalu- 
ating the transport coefficients we will use the newly 
obtained matrix element for gg — > ggg process [||. 



The transport coefficients for QCD matter has 
been evaluated in [Tol - [T3 | . The calculation of the vis- 
cous coefficients within the ambit of diagrammatic 
approach of quantum field theory, along with its lim- 
itation has been discussed in [13| . Recently pQCD 
approaches [l44l8| have been used to calculate rj/s. 
Evaluation of rj/s for a gluonic plasma by Xu and 
Greiner(XG) indicates that the contribution from gg 
— > ggg is 7 times as large as that from gg —> gg. 
This brings the value of rj/s down to the AdS/CFT 
bound ~ 1/4tt [H, 

when the strong coupling con- 
stant, a s = 0.6. ' The GB formula [2(|(see also [U) 
for the gg — » ggg matrix element squared is used 
in [l5[ . However, we have shown recently that at 
the lower temperature domain the GB formula re- 
ceives a significant correction. The ratio of matrix 
element squared with 9] (henceforth will be denoted 
by the subscript DA) and without §(j (henceforth 
denoted by the subscript GB) the correction term is 
given by: 
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where g is the colour charge, N c is the number of the 
colour, mo ~ <?T is the thermal (Debye) mass of the 
gluon, fc^ is the transverse momentum of the emitted 
gluon and q± is the transverse momentum of the 
exchanged gluon. Following our previous work Q 
we depict the magnitude of the correction, R c in 
Fig. [1] The values of q± and s are same as the values 
taken in Ref . @ . It is observed that for large values 
of a s the corrections to the GB matrix element is 
significant. Therefore, it is expected that the values 
of energy loss, rj/s and (/s will also be affected by 
the correction term in the lower temperature (higher 
coupling) domain. 

Before discussing the bulk and shear viscosities we 
estimate the effects of the correction term to the ra- 
diative energy loss mechanism of partons propagat- 
ing through QGP which is measured experimentally 
thorough the nuclear suppression factor [22j in heavy 
ion collision. To evaluate the radiative energy loss 
we start with the soft gluon distribution, which can 
be written as [9| 
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where k — (fco, fcj^, ^3) is the four momenta of the 
emitted gluon, q — (go, 93) is the four momenta 
of the exchanged gluon and Ca = 3 is the Casimir 
invariant of the SU (3) adjoint representation. The 
run in Eq.|3]is the Debye mass required to shield the 
infra-red divergence. We use the above spectrum of 
the soft gluon to evaluate the radiative energy loss, 
dE/dx of gluons. The minimum value of the momen- 
tum of the exchanged gluon in the process: gg — > ggg 
sets time (length) scale for the formation time of the 
emitted gluon. If the formation time is comparable 
to or larger than the mean free time (path) then the 
scattering of the gluons from the successive scatter- 
ers in the medium can not be treated as independent 
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FIG. 1: The variation of the quantity, R c (Eq. [l} with 
coupling constant. 



and the scattering amplitudes between two adjacent 
interaction may interfere destructively leading to the 
the suppression of the emission process - called LPM 
suppression. The LPM suppression has been taken 
into account by including a formation time restric- 
tion on the phase space of the emitted gluon in which 
the formation time, Tp must be smaller than the in- 
teraction time, r = T^ 1 being the interaction 
rate. The radiative energy loss of heavy quark can 
be given by [23l |: 



dE 
dx 



Ye = t ^ .6 



(5) 



where e, the average energy per collision is given by 

dn„ 



e = <n q k > 



drjd fcj 



drjdrk 
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where tf = coshr}/k_i_. The dE/dx is evaluated with 
temperature dependent a s . The variation of a s (and 
hence m 2 D ~ a s (T)T 2 ) with T has been taken from 
Ref. [24| . In Fig. [2] the variation of radiative en- 
ergy loss with T has been depicted for the process 
gg — > ggg. The solid line (dotted line) represents the 
energy loss when DA (GB) gluon multiplicity distri- 
butions are used. To emphasize the importance of 
the corrections to GB formula we display the ratio, 
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in the inset of Fig. [5] It is observed that the cor- 
rection to the gluon spectrum, which leads to the 
energy loss is appreciable for lower temperature do- 
main. This may affect the suppression of high px 
partons in QGP and the elliptic flow of the matter 
formed at RHIC and LHC [25( energies. 
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FIG. 2: The variation of energy loss with temperature 
for the process: gg —¥ ggg. Inset: The Variation of Rel 
(Eq. [7j) with temperature. 



We calculate the quantity, rj/s for a gluonic system 
for the pQCD processes: gg —> gg, gg — > ggg and 
999 ~^ 99- The r\ is evaluated using the procedure 
outlined in [ljj (see also [26| for details): 
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where Bij(k) = B(k)(p 1 k : ' — is defined through 

the infinitesimal deviation from the equilibrium 
value of the gluon phase space density. For the 
process: gg — > ggg we use the matrix element ob- 
tained recently in [si[. The variation of r)/s with 
s — 16 x 27T 2 T 3 /45 as a function of a s is depicted in 
Fig. [3J The ?y/s is quite large at low a s because for 
weakly interacting system the momentum transfer 
between the constituents become strenuous which 
give rise to large rj. However, with the increase in 
the coupling strength the momentum transfer gets 
easier as a result the shear viscosity reduces. The 
results indicate that for large a s the quantity, rj/s 



approaches the AdS/CFT limit. However, in such 
a scenario the non-perturbative effects may become 
important. This can be verified by performing a lat- 
tice QCD based calculations (which include the non- 
perturbative effects) for pure SU(3) gauge theory. In 
fact, such calculation of r\j s has been done in (27j 
and it is found that the value is close to AdS/CFT 
bound. 




FIG. 3: Variation of rj/s with strong coupling constant. 



As mentioned before the bulk viscosity, which is 
connected with the trace of the energy momentum 
tensor through Kubo's formula, will be non-zero for 
a system where the conformal symmetry is broken. 
Lattice QCD calculations indicates non-zero £ for 
a gluonic plasma [5( due to purely quantum effects 
(trace anomaly) (see also [H, |29| for QGP and [3(3] 
for pions) for temperatures around T c . Physically, 
the bulk viscosity appears in the processes which are 
accompanied by a change in the volume(i.e. in den- 
sity) of the fluid. In compression or expansion, as in 
any rapid change of state, the fluid ceases to be in 
thermodynamic equilibrium, and internal processes 
are set up in it which tend to restore the equilib- 
rium. But the processes which drives the system 
toward equilibrium are irreversible associated with 
the increase in entropy and therefore involve energy 
dissipation. Hence, if the relaxation time of these 
processes is long, a considerable dissipation of energy 
occurs when the fluid is compressed or expanded and 
this dissipation must be determined by the bulk vis- 
cosity [31(. We evaluate the bulk viscous coefficient 
with the following formula (see Refs. |32h36| ): 



deg 

rp 
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where deg is the statistical degeneracy for the glu- 
ons, f p is the Bose-Einstein distribution for the glu- 
ons, r is the interaction rate, evaluated using the 
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techniques similar to the one outlined in [17j | with 
the matrix elements Mgg^g g and Mgg^ggg for the 
processes gg — > gg and gg — > ggg respectively, c 2 is 
the velocity of sound and 5 c 2 = (1/3 — c^). The value 
of the velocity of sound, c s for a massless system in 
equilibrium is 1/ -\/3, therefore, the results indicate 
that the bulk viscosity vanishes for a massless system 
in equilibrium. As 5c 2 is a measure of the deviation 
from conformal symmetry (for massless system) the 
£ increases with 5c 2 . 

We evaluate the bulk viscosity with the temper- 
ature dependent c 2 [13] • In Fig- IH the ratio, Q/s 
is depicted as a function of strong coupling. The 
results shown here contain temperature dependence 
thermal gluonic mass. The most striking observa- 
tion one can make here is the completely different 
kind of variation of rj/s and Q/ s with ct s . While CJs 
increases with a s [18j, the rj/s reduces with it jl4| . 
At small a s the bulk viscosity is negligibly small. 
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FIG. 4: The variation of f/s with the strong coupling 
constant. The temperature dependence gluon mass is 
included here. 



We have evaluated the shear and bulk viscosities 
for a gluonic system including the pQCD processes: 
gg — > gg and gg — > ggg. The matrix element for the 
later processes is taken from We find that the 
value of rj/s approaches the AdS/CFT lower bound 
for large a s . The value of rj/s = 0.12 obtained here 
at a s = 0.3 is within the limit extracted from the 
analysis of elliptic flow of matter formed in nuclear 
collisions at RHIC energy (38|. The value of £/s ~ 
0.15 for a s ~ 0.3, which is comparable to rj/s at the 
same value of a s . 

For a s = 0.47 where the corrections to GB for- 
mula is large, we get the AdS/CFT lower bound 
rj/s = 1/47T. Now if we want to understand what is 
the magnitude of rj/s realized in the partonic mat- 
ter expected to be formed in heavy ion collisions at 
RHIC and LHC energies, then one possible way is to 
first estimate the values of temperatures attained in 
these collisions and then needs to know what is the 
value of the a s corresponding to these temperatures. 
The variation of strong coupling with temperature 
may be taken, for example, from (24|, for this pur- 
pose. The typical values of temperatures which can 
be achieved in heavy ion collisions at RHIC and LHC 
energies are ~ 300 MeV and 700 MeV respectively. 
The magnitudes of a s at these temperatures are of 
the order of 0.23 and 0.17 respectively. From the 
variation of rj/s with a s (Fig. [3]) we conclude that 
the typical values of rj/s which may realized at RHIC 
and LHC collisions are 0.177 and 0.29 respectively, 
which is above the AdS/CFT bound but close to 
the value obtained from the analysis of experimen- 
tal data at RHIC m. 
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